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We show that due to entanglement, quantum fluctuations may differ significantly 
■ from statistical fluctuations. We calculate quantum fluctuations of the particle num- 
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briefly discuss the possibility of measuring them. We find that energy fluctuations 
have a non-extensive nature. 

I. INTRODUCTION AND SUMMARY 

In 1964 Bell showed, using his famous inequality, that a system prepared in an entangled 
state can not be described by local, classical physics [1]. Since then, different kinds of 'Bell- 
type-inequalities' have been found and studied. The special nature of an entangled state has 
lead to many developments in quantum teleportation, quantum games, and even in quantum 



black bole physics Q. In this pa ><t wc show thai entanglement is also relevant to the nature 



of quantum fluctuations (See, 



4j for interesting recent discussions of the subject.} 



We discuss a non-relativistic many particle system in a volume Q. We divide the whole 



^ ■ volume f2, into some sub-volume V, and its complement V. The system is initially prepared 
in a pure state (x\ip) defined in Q. We consider states which are entangled with respect 
to the Hilbert spaces of V and V. This means that \ip) can not be brought into a product 
form = ® |^) 2 in terms of a pure state \ip)i that belongs to the Hilbert space of V, 
and another pure state \ip)2 that belongs to the Hilbert space of V. Consider an operator 
O v that is restricted to the sub- volume V. We shall show that the quantum fluctuations of 
O v : (A0 V ) 2 = ((0 V Y) — {(O v )) 2 , behave differently from standard statistical fluctuations 
due to the entangled nature of the initial pure state 

Specifically, we consider energy fluctuations in a sub-volume of a system of cold bosons 
of mass m in its ground state , i.e., a Bose-Einstein condensate (BEC). The BEC is assumed 
to be in a periodic box of size L d (d being the dimension of space.) Had the fluctuations 
been thermal, we would have expected them to be extensive, and so, to scale as the volume: 
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TY&ce(p(AE v ) 2 ) oc V. Instead, we find a different behaviour, for example, the energy 
fluctuations in a sub- volume of size a d scale as the surface area of the sub-volume 

Ag2 = 4^w*y 1+ ,., 

m 2 L 3+d 

Here n 2 Q 2 is a high momentum (UV) cutoff that needs to be introduced in order to render 
the energy fluctuations finite. The UV cutoff is determined by the effective width of the 
boundary that separates the sub-volume from the rest of the system. The dots represent 
sub-leading terms in an expansion in 1/Q. We shall discuss this cutoff shortly. 

The area-scaling is not a special feature of the energy fluctuations in the ground state. 
General arguments [5] show that it occurs for a large class of operators whose two-point 
function is short ranged. It has also been shown that the area scaling is robust to changes 
in the cutoff scheme and does not depend sensitively on the nature of the pure state. 

The energy fluctuations exhibit a power law dependence on the high momentum cutoff. 
This may seem strange and counter intuitive. One may ask why should a macroscopic 
quantity depend so strongly on the details of the physics that are relevant only to very 
short distances? Our explanation for these power law divergences relies on the uncertainty 
principle [l^; a measurement of the momentum fluctuations of a localized particle is formally 
divergent. Smearing the position of the particle will render the momentum fluctuations finite 
and inversely dependent on the smearing scale. The divergence in the energy fluctuations 
is of similar origin: measuring the energy fluctuations in a volume with a sharp boundary 
yields an infinite result. Introducing a high momentum cutoff is equivalent to smearing the 
boundary of the sub- volume. The width of the boundary is inversely dependent on the cutoff 
scale. 

A measurement of energy fluctuations in the laboratory may shed some light on these 
issues. Experiments in BEC type systems seem to offer a unique opportunity for measuring 
entanglement induced fluctuations. To perform such an experiment one needs to be able 
to make repetitive measurements of a particular sub- volume of the condensate, and prepare 
the condensate in the same initial state for each measurement. Alternately, one may make 
measurements of several sub- volumes in the same condensate using tomographic methods. A 
verification of the explicit dependance on the high momentum cutoff of energy fluctuations 
is not only interesting in its own right, but may also have implications concerning Unruh 
radiation and black hole physics [6j. In the BEC context one may guess that the UV 
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momentum cutoff is inversely proportional to the resolution of the measuring device. If 
so, our result in eq.flTJ implies that a better detector resolution will lead to larger energy- 
fluctuations. 

Much has been said about the area law of entanglement entropy (see for example M, M) 

nnnn 

and black holes [9|, HO, 1111, 1121]. To understand the possible relationship of entanglement 
induced fluctuations to black hole physics it is useful to consider as an example the vacuum 
state of some field(s) in the space-time of an eternal black hole. An observer outside the 
black hole will make measurements that are restricted to the exterior of the event horizon. 
Let us take for example the measurement of energy fluctuations by this observer. They are 
divergent, depend on a power of the high momentum cutoff, as is the entanglement entropy, 
and scale as the surface area of the horizon. The divergence originates from distance scales 
which are very close to the horizon. Technically they arise from the divergence of the two- 
point function of the energy density at short distance scales. This is the same mechanism 
which generates the divergences in the case of cold bosons. The understanding that we 
have gained on the nature of the divergences in the non-relativistic setting provides a new 
perspective on the nature of the same divergences that appear in the calculation of energy 
fluctuations in black hole spacetimes. 

In addition to the energy fluctuations, one may also consider fluctuations in the particle 
number (see [13J] for a recent discussion.) We give an analytic expression for the fluctuations 
of the particle number in various potentials in section IIHI The particle number fluctuations 
are different from the energy fluctuations in that they do not depend on the high momentum 
cutoff. This is because the two point function of the number density operator does not 
diverge strongly enough at short distances 5|. Yet, we show that the quantum correlations 
in the initial state generate fluctuations which are different from the ones that are expected 
in a classical setup. 



II. GENERAL SETUP 



We shall describe the general method of evaluating fluctuations of operators restricted to 
subvolumes in the case that the whole system is in a pure state We will consider the 
case that the pure state has a fixed number of particles and fixed momentum. We discuss 
systems whose Hamiltonian is time-independent and have some fixed boundary conditions 
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(BC), for example, periodic or Dirichlet. For the general discussion we do not specify the 
BC, whose significance will be discussed later when we consider in detail fluctuations of two 
specific operators: the particle number operator, and the energy operator. 

We start with a general operator density O(x). It can be expressed in terms of the 
particle field ijj(x,t) = ^ i aiUi{x)e~ lEit as follows, O(x) = ^(x,t)Oijj(x,t). The functions 
basis of eigenfuctions of the Hamiltonian. For example, the energy operator is given 
by E = — ^;V 2 + V, and the particle number operator by N — I. 

We consider a pure state \ip) with n particles of momentum k. Evaluating the operator 
O v = J v 0(x)dV in this state gives 

<^|o» = m J J2^v u i(x)du p (x)dv\^} 

i, P 

= J ul(x)du k (x)dV = nO Kk . (2) 

The fluctuations of the operator O v , (^|(AO v ) 2 |^> = (i/j\(O v ) 2 \i/j) - ((iJ>\O v \ip)) 2 , are given 
by 



^kk 



= n^2d ks d sk . (3) 

s^k 

One may generalize this result to a system with T particle modes, such that there are 

(a) ) n * 

ni {1 <i <T) particles with quantum numbers k iy \ip) = Hj =1 —^=-\0). The most general 
state will be a linear superposition of such multi particle states. In this case we have 

T 

^\O v \^) = Y,n l d ktM , (4) 
i=i 

and 

(° V ) V) = Yl KidkisOski + n i n eOk lke O keki . (5) 

In order to study the properties of the fluctuations we need to evaluate the matrices 
O ks = J u k (x)Ou s (x)dV . We shall consider two special cases: the number of particles 
operator and the energy operator. For the number operator we have N — I. Then 

Y^N ks N sl = V / u* k (x)u s (x)dV x [ u*MMv)dVy 

= [ ui(x)6(x-y) Ul (y)dV x dV y = N kl . (6) 
Jv 
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So for one particle mode we get 

(AN V ) 2 = nN kk - nN 2 kk = (N v )(l - N kk ). (7) 

For the energy operator, we have: 

E ks = [ ut(x)Hu s (x)dV 
Jv 

= E s N ks (8) 
So that for the case of one particle mode we get: 

(AE v ) 2 = nY,E k E s N ks N sk . (9) 

s^k 

III. RELATIVE PARTICLE NUMBER FLUCTUATIONS 

We proceed to consider in detail the relative particle number fluctuations in the case that 
the initial state contains a single particle mode: (AN V ) 2 / (N v ) = 1 — N kk . 

A. Free particles in a box 

The simplest case to consider is a one dimensional box of size L, with periodic BC. The 
solutions to the Schrodinger equation are of the form u n {x) = p^e 2 ™ 2 ^, with n E Z. One 
obtains N n ^ n = We also give the result for the case m ^ n for future reference, 

7V n , m = 1/L f e~ 2mnx/L e 2mmx,L dx (10) 
Jo 

= — e Mm-n)a/L ^faOm - Tl)a/L). 

7r(m — n) 

One then has (AN V ) 2 /(N V ) = 1 — a/L. Generalizing to d-dimensions, we get 

(AN v ) 2 /(Ny) = 1 - {a/L) d . (11) 
The eigenstates for free particles in a box with Dirichlet BC (Dirichlet box) are: 

u n (x) = a/2/ L s'm(nnx / L) , 

and so 

-~i 1 /sin((m — n) TV a/L) sin((m + n)na/L)\ 

Nn,m = ~ I : I • { 12 J 

7T V m — n m + n I 
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FIG. 1: Shown is a graph of the relative particle number fluctuations (ANv) 2 /(Ny) in a sub- 
volume of size a of a one dimensional box of size L with Dirichlet boundary conditions versus the 
relative size of the sub-volume a/L. Each color corresponds to a different initial state as specified 
by the legend. 



The diagonal matrix elements are 

at a sin (2nkf ) 
Nk > k ~L 

In d- dimensions one has 

^^.(..tM) 4 . (13) 

The first five levels for a one dimensional system are shown in Figure [D The results for the 
Dirichlet box are qualitatively similar to the results of the periodic box, the main difference 
being the oscillations that are superimposed on the factors of a/L. These oscillations are 
purely quantum and do not appear in an equivalent microcanonical ensemble as we show 
next. 

We would like to compare our results in equations (fTTl) and ffT3"|) to their classical analogs. 
We show that the case of the periodic box gives identical results to the case of a classical 
uniform distribution of particles in a box. The case of Dirichlet box is different from the 
classical example and exhibits interesting quantum oscillations. 

Let us consider a classical system with n particles, uniformly distributed in a box of size 
L and calculate (AN V ) 2 and (N v ). The probability for finding a particle in a given position 
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is : P{x)dx = 1/Ldx. Therefore the probability of finding a particle in a region of size 
a < L is p = a/L. The probability of not finding it there is g = 1 — a/L. If there are N 
particles and they are distributed so that each one has a uniform probability to be found in 
a given position, then the probability of finding exactly N a particles inside a region of size 
a is binomial 

P(N a \N) = ( N | p"V-<V (14) 

Therefore the mean number of particles is (N a ) = Np, and the variance is AiVf = Npq. So 
that 

^-4-1-- (15) 

This result generalizes in an obvious way to d-dimensions = 1 — (f ) d , which is exactly 
the result in eq. ljTTj) . The behavior exhibited in eq. liTBT) is quite different from the one 
exhibited in eq. (lT3i) for the Dirichlet BC case where the correlations may induce a sinusoidal 
variation in the fluctuations (see Figured!) 



B. The harmonic oscillator 

The BEC condensate may be described as a collection of bosonic particles in a harmonic 
oscillator potential (see [14j for a review), so it is interesting to consider the particle number 
fluctuations in such a setup. The eigenstates for a one dimensional harmonic oscillator are 
given by 



«„(•'■) 4= 2 " 2 ( 4=] ' /U '/•'')' 



'71! W 71 ", 

where we have defined r\ = ^Jmuo/h. Therefore, if our sub-volume is the region a < x < b, 
then for m/nwe have: 



~ y/n\m\ sr-^ 2' / _„ 2 „2 TT . . _„2 fc 2 TT . ,.\ 



and for m = n we have: 



I r — — ^ ^ ^22 2 l2 



2"v^t^«(«-0Kn-0! 



=o 

1 



+ -(erf (776) — erf(r^a)). 
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FIG. 2: Shown is a graph of the relative particle number fluctuations (AiVy) 2 /(iVy) for the case 
of a harmonic oscillator of mass m and frequency uj. The sub- volume is in the region —a<x<a 
and the fluctuations are plotted against the relative size of the sub-volume in units of WrruJ/H. 
Each color corresponds to a different initial state as specified by the legend. 

For the case that the sub- volume is symmetric about the minimum of the potential —a < 
x < a the expressions simplify and we find that for the lowest energy state ((ANy) 2 ) / (Ny) = 
1 — erf (rja). This, and a plot of the first four excited modes (in one dimension) is shown 
in Figure [2J For the lowest level of a d-dimensional symmetric harmonic oscillator we have 
(AAV)V(AV) = 1 - (erf(r/a)) d . 

The particle number fluctuations in the lowest energy state of a harmonic oscillator have 
some similarities to the fluctuations of free particles in a box. The "size" of the whole 
system, L in this case, is approximately L ~ 1/rj (recall that r] = \J mw / K) , and for rja < 1, 
erf(r/a) ~ 2/yGra/L. It follows that {AN V ) 2 /(N V ) ~ 1 - (a/L) d . However, in the higher 
energy states, the "size" of the system increases since the wave function spreads more and 
more as the energy becomes higher, and the particle number fluctuations exhibit a different 
behaviour. 

IV. ENERGY FLUCTUATIONS 

The derivation of the result for the energy fluctuations from eq. (Q is not as tractable as 
the derivation of the particle number fluctuations. We evaluate the energy fluctuations in 
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a cubic sub-volume of a d-dimensional periodic box and show that the energy fluctuations 
scale as the surface area of this sub-volume. It should be clear that result is quite general, 
and that similar results should be expected for other setups. There is also ample additional 
evidence that this scaling behavior is quite robust 

The expression for energy fluctuations is given in eq.([9]). We start by evaluating 
J2 s E k E s N ks N sk . Using eq.flTD), we get: 

£ = ^ t W ->i)- ("J 

s — oo 

Igl 

This expression is divergent, so we introduce an exponential UV cutoff on the sum e ^Q' 2 . 

As shown in the Appendix, the cutoff introduces a smoothing of the boundary of the 
sub-volume. So far we have considered a perfectly sharp boundary < x < a and did not 
encounter any problems. Equivalently, we have defined the sub-volume by integrating with 

{1, < x < a a L 
, so J dx = J dxH(x — a). 
0, x>a 

The smoothing of the boundary is over a region x ~ L/Q, and equivalently the smooth 
step function changes its value from 1 to over a range ~ L/Q. The specific cutoff scheme 
determines the detailed behaviour of the smoothed step function over the region where 
it changes significantly. However, all the cutoff schemes induce very similar qualitative 
smoothing. The advantage of using exponential cutoff is that the expression for the energy 
fluctuations can be evaluated analytically with the cutoff. 
With the cutoff eq. (fT6l) becomes 

E E * E AsN sk -"^E 7^ sin V( fc " ( 17 ) 

s — oo ^ ' 

This series may be evaluated analytically. We find 

OO 9 

v — v S I s 

n 2 dE(Q, k,x) = \ — — — —sm 2 (ir(k-s)x)e~^ T Q 2 ~ = 

\K S) 

s=—oo v ' 

k coth ( 2„ 2 ) sin 2 (to) *• 1 1 1 i 

+e *■ « -(2Li 2 (e " « J— Li 2 (e ^ « )— Li 2 (e « )), 



cos(27rx) - cosh ^^2- j 2 



where x = a/L, and Li 2 is the dilogarithm function. 
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FIG. 3: Logarithmic plot of energy fluctuations versus a/L for a one dimensional box with periodic 
boundary conditions. Here, the UV cutoff tt 2 Q 2 is 10000, the state is the 20th excited state (k = 20), 

E 2 

and the mass has been chosen such that = 1. 

Therefore, the energy fluctuations for n particles in one dimension are given by 

E 2 

AE 2 = n-jdE(Q, k, a/L). (19) 

A plot of AE 2 is shown in Figure [31 We wish to emphasize several points: first, at x = 
(a = 0), and x — 1 (a = L), the energy fluctuations vanish. This is as expected for a volume 
of zero size, and the volume of the whole box respectively. However, these two points are 
discontinuous in the Q — > oo limit. In the vicinity of these points dE(x) rises sharply and 
then remains at an almost constant value. This can be seen by expanding in powers of 1/Q 
for < x < 1. We find 

dE(Q, k, x) = Q 2 + \k 2 {\ - (1 - 2x) 2 ) - A + 0(_L). 

If we consider fluctuations for a low energy state (though not the ground state), we have 

AE 2 = nEg + Oil). (20) 

Where, again, the expansion parameter is 1/Q. Note that the UV divergence of the energy 
fluctuations is linear in the cutoff parameter. Our choice of n 2 Q 2 as the cutoff was made to 
bring eq. (l20~!) into a convenient form. 
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In the d- dimensional case, we have: 

d 



n i 1 v 



3 27r(m; —rii)j- 



l sin 7T m,- — n,- 



while 



£ 2 47r 2 s 2 
2mL 2 



Now, 



E *wu*» = E ■ ■ ■ E E A - ■ s,)a/i)s? 



' 2mL 2 7i 

s si sn j i=l 

Defining the quantity M(si) = — ' " ! " - '" — 



;2i) 



^EEEni^w^*-^ 



( fc ._ s .)2 > we can rearrange the sum in eq. (l2"T]) as follows, 

d-l~ 



si 



*2 



k=l s k j=l 

Now, according to eq. (fT8|) . Y^ si M(si)sf = n 2 dE(Q, k,x) and Y^ S2 M(s 2 ) is the one dimen- 
sional N mm operator which was shown to be equal to a/L. Therefore, the energy fluctuations 
are given by 



AE = ndE, 



h 2 4:Tr 2 dE(Q,k,x) fa 
' 2mL 2 



d-l 



(22) 



If Q is large and a/L is small this reduces to 

AE 2 = nE k E Q d + 



which has the claimed area scaling behavior 16]. 



V. CONCLUSIONS AND DISCUSSION 



In this paper we calculated fluctuations of the particle number operator and the energy 
operator restricted to sub- volumes. The whole system was in a pure state with fixed particle 
number and fixed momentum. We discussed systems whose Hamiltonian is time-independent 
with periodic or Dirichlet boundary conditions. In general the fluctuations differ from their 
classical counterparts. 

Even for a single non-relativistic particle in a box, the energy fluctuations in a sub- volume 
with sharp boundaries are divergent and scale with the area of the boundary of the sub- 
volume. This behaviour depends weakly on the specific state of the system, its boundary 
conditions, or on the specific Hamiltonian. 
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In addition to scaling with the area of the boundary the energy fluctuations diverge with 
the UV cutoff. Technically this cutoff appears when we sum over energy eigenmodes: by 
introducing a high-momentum cutoff on the Fourier modes in the sum, the energy fluctu- 
ations become finite and proportional to the value of the UV cutoff. We have shown that 
introducing the cutoff is equivalent to smoothing the boundary of the sub-volume such that 
its width is inversely proportional to the cutoff. 

One might have thought that the fact that the divergences are associated with the UV 
modes suggests that the origin of the divergence requires some additional information on the 
UV behaviour of the system. We believe that the resolution lies elsewhere: it seems that the 
divergence of the energy fluctuations is a result of our assumption regarding sharp boundary 
conditions of the subvolume. Such sharp boundaries would be impossible to construct in 
practice. This is similar to an attempt to compute momentum fluctuations of a particle 
which is sharply localized: if the position of the particle is known with infinite precision 
then the momentum fluctuations formally diverge. Smearing the position of the particle will 
render the momentum fluctuations finite and inversely proportional to the smearing scale. 
It seems as if the divergence in the energy fluctuations is of similar origin, i.e., measuring 
the energy fluctuations in a volume with an infinitely sharp boundary yields an infinite 
result. Introducing a high momentum cutoff is equivalent to smearing the boundary of the 
sub-volume. Put differently, the width of the boundary is inversely dependent on the cutoff 
scale. A similar effect appears in relativistic field theories js| and it would be interesting 
to strengthen these suspected ties with the uncertainty principle by performing additional 
calculations. 

For the case of particle number fluctuations we found that the differences between the 
quantum and classical cases are more subtle. This is related to the fact that number fluc- 
tuations do not formally diverge even for an infinitely sharp boundary of the sub-volume. 
Nevertheless there are differences between the classical fluctuations and the quantum one 
which depend more strongly on the boundary conditions, the nature of the global quantum 
state and the Hamiltonian of the system. 
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APPENDIX A: FOURIER TRANSFORMS OF SMOOTHED STEP FUNCTIONS 

The basis of eigenfunctions in a one dimensional periodic box of size L = 1 is ip n = e 2mxn 
where n is an integer — oo < n < oo. Any periodic function can be expanded in this basis, 



and the coefficients are given by 



fn 



/ dxf(x)e 
Jo 



-2-Kixn 



(Al) 



(A2) 



The (periodic) step function is defined as H(x — a) 
lowing expansion coefficients 



1 < x < a 



dxH(x — a)e 



-2mxn 



x > a 



sm[7ranj 
irn 



and has the fol- 



(A3) 



The step function has a useful analytic representation, 

g— 2ni(a— x) ' 



H(x — a) 



a + — log 

2m 



I g+2-7ri(a— x) 



+ log 



1-e 



-2nix 



J £,+2nix 



(A4) 



Defining z = 1 — e 2?n(a x > and 



w 



-2nix 



we observe that logarithms depend only on 



the ratios z/z and w/w, hence they depend only on argz and argw. The steps occur when 
the branches of the logarithms change abruptly at x = and x = a when the absolute value 
of w or z vanish. 

a 

The integral J f(x)dx can be treated as a convolution with a step function and hence 



J f(x)dx = J dxH(x — a)f(x) =^^H*f n . (A5) 

n 

In eq.( |T6i) appears a sum of the form H*_ k H s _kE s Ek which results from a product of two 

s 

step functions and two energy operators. 

Introducing a UV cutoff in the sums that appear in eq. (1A5l) or in eq. (fT6|) amounts to 
suppressing the coefficients H n beyond a certain large n c , H n — > 0, n > n c . An exponential 
cutoff amounts to replacing H n — > H n e~^ n ^ nc . To see explicitly the effect of an exponential 
cutoff we may evaluate the position space representation of the cutoff expansion 



H(x — a) = H n e «■= = a + 



2m 



log 



1 — e n 



1 2m(a— x) 



-i — — +2-7ri(a— x) 
1 — e "c v ' 



+ log 



1 — - — 2irix 
-i - — +2irix 

1 — e ™c 



(A6) 
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The branches of the logarithms now change over a distance Ax ~ 1/ ^fn~ c near the points 
x = and x — a. For an interval of size L rather than 1 the change is over a range 
Ax ~ Lj Tjn~ c . 

We have verified that other cutoff schemes yield qualitatively similar results by performing 
numerically similar calculations for different cutoff schemes and for different BC. 
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